201711.02295v1 


chinaXiv 
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collected with the Beijing Spectrometer (BESIII). For the first time, the branching fractions of xc > 


K^ K, xco — a; (1320)? 


F /a5(1320)x? and xc — p(770)** are measured. Here K*K denotes both 


K** K* and its isospin-conjugated process K*°K° + c.c., and K* denotes the resonances K* (892), 
K3(1430) and K3(1780). The observations indicate a strong violation of the helicity selection rule 


in Xe2 decays into vector and pseudoscalar meson pairs. 


The measured branching fractions of 


Xe2 > K* (892) K are more than 10 times larger than the upper limit of xc > o(770)^ 7, which 
is so far the first direct observation of a significant U-spin symmetry breaking effect in charmonium 


decays. 


PACS numbers: 13.25.Gv, 12.38.Qk 


The helicity selection rule (HSR) is one of the 
most important consequences of perturbative quantum 
chromodynamics (pQCD) at leading twist accuracy. In 
the charmonium energy region, although there are ob- 
servations that pQCD plays a dominant role, there are 
also many hints that non-perturbative mechanisms can 
become important . Exclusive decays of the P-wave 
charmonium state Xe2 — V P, where V and P denote 
light vector and pseudoscalar mesons, respectively, are 
ideal for testing the HSR and pinning down the mecha- 
nisms that may violate the leading pQCD approximation. 


Another reason the decays of x.9 — V P are of great 
interest is that this process is ideal for probing the long- 
range interactions arising from intermediate D-meson 
loop transitions. As pointed out in Ref. i, if the in- 
termediate D-meson loops provide the non-perturbative 
mechanism to violate the HSR, this can be identified 
by the measurements of y, — K*(892)K and ye — 

p(T10) rF. 

In this Letter, we present a partial wave analysis 
(PWA) of the process xc > KK7 (denotes K* K ^x? 
and KsK-7-*) and a measurement of xc — «^s m. 
We have two 1(3686) samples of (106.8 + 0.8) x 10° 
(160 pb~!) [8| and (341.1 + 2.1) x 109 (510 pb-!) [gj 
events collected in 2009 and 2012 by BESIII , respec- 
tively. Only the 2009 data sample is used in the anal- 
ysis of y,» ini KK, and the full data sample is used 
in Xe9 > m ^7 m? since it has a smaller branching frac- 
tion. An ed sample of about 44 pb^! taken at 
Vs = 3.65 GeV is utilized to investigate the potential 
background from the continuum process. A sample of 
Monte Carlo (MC) simulated events of generic (3686) 
decays (inclusive MC sample) is used to study back- 
grounds. The optimization of the event selection and the 
estimation of physics backgrounds are performed with 
Monte Carlo simulations of (3686) inclusive/exclusive 
decays. 

The xez candidates, produced in (3686) radiative oy 
cays, are si i Shu from the final states K+K~7° 
KgK-m«*, and ntr’. Each charged track is re- 
quired to have a polar angle 0 in the main drift cham- 
ber (MDC) that satisfies | cos0| < 0.93, and have the 
point of closest approach to the ete~ interaction point 
within 10 cm in the beam direction (|V;|]), and 1 cm 
in the plane perpendicular to the beam direction (V,). 


The energy loss dE/dx in the MDC and the information 
from the time-of-flight (TOF) system are combined to 
form particle identification (PID) confidence levels (C.L.) 
for the 7, K, and p hypotheses, and each track is as- 
signed with the hypothesis corresponding to the highest 
C.L. The Kg candidates are reconstructed from two op- 
positely charged tracks with loose vertex requirements 
(\V.| < 30 cm and V. < 10 cm) and without PID (as- 
sumed to be pions). Then the candidate with invari- 
ant mass closest to the Ks nominal mass and the de- 
cay length provided by a secondary vertex fit algorithm 
greater than 0.25 cm, is selected for further study in the 
decay (3686) > yKsK*+x*. The candidate events are 
required to have two charged tracks with zero net charge, 
where the tracks from the Kg candidate are not taken 
into account. Two nuni and one kaon are required for the 
decays 1/(3686) — »yr* 1-1? and (3686) > yKsK=nF, 
respectively, and no PID requirement is applied for the 
decay (3686) — yK*K~nx°. The photon candidates 
are required to have energy lod than 25 (50) MeV in 
the Electromagnetic Calorimeter (EMC) barrel (end cap) 
region |cos0| < 0.8 (0.86 < |cos0| < 0.93), and have an 
angle relative to the nearest charged tracks larger than 
10°. To suppress electronic noise and energy deposits 
unrelated to the event, the EMC cluster time must be 
within 700 ns from the event start time. At least three 
i one photons are required for the decay (3686) — 
yr*n n1?9/K* Kn? and (3686) 2 4KsK 7, respec- 
tively. 


A fit with four kinematic constraints (4C) enforcing 
four-momentum conservation between the initial (3686) 
and the final state is performed for each process. If there 
are more photons than required in one event, all pos- 
sible combinations of photons are considered and only 
the one with the least xj, of ue kinematic fit is re- 
tained for further analysis. The x26 is required to Da less 
than 80 and 60 for the decay 4/(3686) > aS 1? and 
1,(3686) — yKsK*rF, respectively. The 7° candidate 
is reconstructed from the two selected photons whose in- 
variant mass is closest to the 1? nominal mass, and sat- 
isfies |My. — M, RI e 10 MeV/c?. For the decay mode 
(3686) + yrtr-7°, a 5C kinematic fit is petonined 
with an additional 7° mass constraint, and x25 < 60 
is required. To remove the backgrounds (3686) — 

070 J p (Jf — U017, l= e, u), the invariant mass of 
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K*K-[n*n- is required to be less than 3.0 GeV/c? 
for the decay (3686) — yKtK~7°/yrtn-n®. In 
the decay mode (3686) — -*m«*m--7?, the 2° recoil 
mass is required to be less than 3.0 GeV/c? to sup- 
press the background (3686) > 519J/w, and Myro € 
(0.7, 0.85) oid is feguired to veto the background 
vV(3686) — wrtn (w — »yn?). 

The K id T PARC mass for the decay (3686) — 
yK*K-n? and (3686) + yKsK+n*F are shown in 
Fig. [{a) and (b), respectively. The xci signals ap- 
pear prominently with a small background. From the 
analysis of the 7(3686) inclusive MC sample and the 
continuum data at ys = 3.65 GeV, the main back- 
grounds are from the decays w(3686) —  m?«?J/w 
(J/y — ata-n°/pp), and (3686) > K1(1270)* KF 
(K,(1270)* — K-n?n9/p(710)* Ks). All of these back- 
grounds show a smooth distribution, and do not pro- 
duce a peak around the xy.; mass region. Unbinned 
maximum likelihood fits are performed to the selected 
candidates, where the x.1,2 signals are described with 
the MC simulated shapes convoluted with a Gaussian 
function accounting for the resolution difference between 
data and MC simulation, and the backgrounds are de- 
scribed with a 2"4 order polynomial function. There are 
1215 and 1176 candidate events for ye — KtK 771° 
and ye —  KgK-m«* within the xe» signal region 
Mice — My] € 15 MeV/c?. Non-s? (K* K ^1? mode 
only) and non-y,-2 backgrounds are estimated with the 
events in the sideband regions, which are also used in 
the PWA as described in the following. The numbers of 
background events are estimated to be 240 and 80 for 
Xe > KKT’ and xc > KsK** , respectively. 

In the PWA, the process xe2 + KK is assumed to 
proceed via the quasi two-body decays, i.e. Xe2 — asm 
and K*K followed by a2 > KK and K* + Kr. The am- 
plitudes of the two-body decays are constructed with the 
helicity-covariant method mi For a particle decaying 
into two-body final states, i.e. A(J, m) + B(s, A)C (0, v), 
where spin and helicity are indicated in the parentheses, 
its helicity covariant amplitude F) ,, is: 


Fyw = M Agus(L089|Jó) (so — v|$6)r"Br(r), (1) 
LS 


z 
2J+ 


partial wave with abid angular momentum L and spin 
S (with z-projection ó), r is the relative momentum be- 
tween the two daughter particles in the initial particle 
rest frame, and By is the barrier factor [13]. The con- 
servation of parity is applied in the equation. Recent 
measurements show that the contributions of higher or- 
der magnetic and electric multipoles in the (3686) ra- 
diative transition to y-2 are negligible, and the E1 tran- 
sition is the dominant process Hence, the helic- 
ity amplitudes. are constructed to paiia the El transi- 
tion relation [14] and parity conservation, namely, Fi» = 


where A — , gts is the coupling constant for the 
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FIG. 1. (color online) Invariant mass distribution of (a) 
Kt Ks? and (b) Ks K*n* for the decay (3686) — yK Km, 
and the corresponding Dalitz distributions (c) and (d) for the 
candidates within |Mrx« — My,4| < 15 MeV/c?. The dots 
with error bars are for data, the blue solid curves are the over- 
all fit results, the red dotted curves are the signals, and the 
green shaded areas are the background. 


V2Fi1 = V6F\.0 and Fo o = 0. The corresponding grs 
are taken as complex values. The relative magnitudes 
and phases are determined by an un-binned maximum 
likelihood fit to data with the package MINUIT iş. The 
background contribution to the likelihood value is esti- 
mated with the events in the sideband regions and is sub- 
tracted [1g]. For the PWA method check, a input data 
is generated with inclusion of all states in the baseline 
solution, and coupling constants are fixed to the PWA 
solution. After the detector simulation and selection cri- 
teria, the same PWA fit procedure is performed, the fit 
results are consistent with that of the input data within 
the statistical errors. 


As shown in the Dalitz plots of Fig.[I(c) and (d), clear 
signals for K*(892) and K3(1430) are observed in the 
Km system. The resonances K*(892) and K3(1430) in 
the Kr system as well as the a2(1320) in the KK sys- 
tem, which has a significance larger than 8e in both de- 
cay modes, are included in the baseline solution. For 
consistency, the A3(1780) signal, which is of a signifi- 
cance larger than 5c in the decay ye2 > KsK-m*, but 
only 2e in Xe > K^ K- «9, is also included. A contri- 
bution from the direct Xe2 — KK three-body decay, 
which is parameterized with a non-resonant component 
with spin-parity J? = 2+ in the KK system, is also con- 
sidered. Other possible excited K* states in Kr system 
and the states in the KK system listed in PDG mm 
which have a significance less than 5c, are not included, 
but they are considered as a source of systematic uncer- 
tainty. The coupling constants for the charge-conjugate 


modes are treated to be the same. 

Figures [2] and [3[a)-(c) show the invariant mass dis- 
tribution and the projection of the PWA for the decays 
Xe2 > K* K-97 and yo > KgK m7, respectively. The 
signal yields for the individual processes with a given in- 
termediate state and the corresponding statistical uncer- 
tainties are calculated according to the fit results. The 
resultant branching fractions for the decays xe > K*K 
and xXc2 — a@2(1320)7 are summarized in Table [I] The 
branching fractions for the processes including charged 
K* intermediate states are consistent between the two 
decay modes, and are combined by considering the cor- 
relation of uncertainties between the two modes [18]. 
The K* isospin-conjugate modes are consistent with each 
other within 20, as expected by isospin symmetry. 
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FIG. 2. (color online) Projections of the fit results onto the 
invariant mass of (a) K*° and (b) KK” in the decay 
Xc2 — Kt K~7°, where dots with error bars are for data, 
the blue solid histograms are the overall fit results, the yellow 
shaded histograms are for the background estimated using xc» 
sideband events, and the contributions from different compo- 
nents are indicated in the inset. 


TABLE I. The measured branching fractions (x10~*) for the 
decays Xe2 — (Kr) K and xe3 — (XK)r. The first uncer- 
tainties are statistical, and the second are systematic (Here, 
K*, K3, K3, and a2 refer to K*(892), K5(1430), K3(1780), 
and a2(1320), respectively). 


Mode K Kn” Kek=nt Combined 
K*'*K* 18+£02£02 14£02+£02 15£01402 
KOIR” — 1320.24 02 = 
KEK 18.2+0.841.6 13.6 + 0.8 + 1.4 15.5+064+1.2 
Ki?K? z 13.0 + 1.04: 1.5 = 
Ki*K* 53405409 59411415 5.40.54 0.7 
Ki?K? = 5.9 1.6 X 1.5 = 

abn? 13.5 + 1.6 + 3.2 = = 

az T — 18.42: 3.3 £5.5 — 


The large branching fraction of xc + K3K is a di- 
rect indication of the significant HSR violation effects. 
Note that the helicity amplitude ratios | F2 o|?/| F1 o|?, es- 
timated with the fitted grs (see Table [II], suggest the 
dominance of F} o in the transition amplitudes. The am- 
plitude Fo contributes to the leading HSR violation ef- 
fects and scales as (Agcp)/m-)® due to its asymptotic 
behavior B, [j. In comparison with the HSR conserved 
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FIG. 3. (color online) Projections of the fit results onto th 


invariant mass of (a) K^, (b) Ksz^ and (c) KsK- in 
the decay xc > KsK+*n* as well as (d) «^a? in the decay 
Xc2 — cm m, where dots with error bars are for data, the 
blue solid histograms are the overall fit results, and the con- 


tributions from other components are indicated in the inset. 


channel xc — VV, which scales as (Agcp/m.)^, the ra- 
tio of xc + K3K to VV is expected to be suppressed by 
a factor of (Agcp/m.)? ~ 0.02 with Agcp ~ 0.2 GeV/c? 
and the charm quark mass me ~ 1.5 GeV/c?. How- 
ever, the measured branching fraction of Xe2 — KiK 
appears to be the same order of magnitude as that for 
Xe > VV [19]. which indicates a significant violation of 
HSR in xc > K3K. 


TABLE II. The measured ratios of helicity amplitude squared 

[F2.0]?/|F1,0|?, where the uncertainties are statistical only. 
KK n” KsK=n* 

Charged K* Charged K* Neutral K* 

0.046 + 0.001 0.042 + 0.019 0.031 + 0.018 


K (1430) 


In the analysis of the decay (3686) — yt 7°, the 
Xc2 signal is extracted by the requirement |M,+,-,0 — 
M4] € 15 MeV/c?. The potential background from di- 
rect e* e^ annihilation is found to be negligible by study- 
ing the continuum data taken at ys = 3.65 GeV. The 
backgrounds from 4(3686) decay are investigated with 
the 1/(3686) inclusive MC sample; the only surviving xc» 
events are those that directly decay to t+ a7 7° without 
any intermediate state. There are also non-%e2 back- 
grounds, which can be estimated by the events in the x. 
sideband regions. Figure [B{d) shows the invariant mass 
of 1*7? for the selected candidates, together with the 
binned likelihood fit results. Here the fit components in- 
clude the p(770)* signal, the direct decay xoa — m^ m 7° 
and the non-x,2 background. The p(770)* signal and 
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the direct Xe2 three-body decay are modeled with the 
MC simulated shapes convoluted with a Gaussian func- 
tion with free parameters. The resonant parameters 
of the p(770)* are set to the values in the PDG [13. 
The fitted signal yields are 14.7 + 8.9 and 63.6 + 13.0, 
and the corresponding resultant branching fractions are 
(0.64 3- 0.39 3- 0.07) x 107? and (2.12: 0.4-£ 0.2) x 107? for 
Xc2 — p(770)*n* and the direct decay xc3 > ata m9, 
respectively, where the first uncertainties are statistical, 
and the second are systematic. Since the statistical sig- 
nificance for the xc — p(170)-77 is only 2.80, the upper 
limit at the 9096 C.L. for the branching fraction is set to 
1.1 x 107? by the method of Feldman-Cousins approach 
with the systematic uncertainties consideration bd 

The uncertainties from the branching fractions of 
v(3686) > 7x02, Ks > ata, no — yy, ant) > 
KK, and K* — K are quoted from the PDG [17]. The 
uncertainty on the number of (3686) events is abont 
0.8% is, . The uncertainties associated with the track- 
ing and PID are 1% for every charged track [21]. The 
uncertainty related with EMC shower reconstruction ef- 
ficiency is 1% per shower [21]. The uncertainties asso- 
ciated with the kinematic fit are estimated to be 0.596 
and 0.696 for the 4C and 5C fit, respectively, by us- 
ing a method to correct the charged-track helix param- 
eters [22]. The uncertainty associated with the Kg re- 
construction is estimated to be 2.596 kg. The uncer- 
tainties related with the 7° selection, the requirements 
on the 7° recoil mass and the w background veto (in 
Xc2 > p(770)*x* mode only) are negligible. 

In the decay xc» — p(770)*7*, the uncertainties due 
to the bin size and the fit range in the fit are estimated by 
repeating the fit with alternative bin sizes and fit ranges. 
The uncertainty due to the shape of Xe2 — ntron’ is 
estimated by replacing the MC simulated line shape with 
a 3'7 polynomial function. The uncertainty due to the 
shape of the background is estimated by changing the 
Xe2 sideband regions. 

In the decay x. + K Kr, the uncertainties due to the 
contribution from K*(1410) and K*(1680) are estimated 
by including these states in the fit. The uncertainties as- 
sociated with the backgrounds are determined by chang- 
ing the y-2 and 7° sideband regions. The spin density 
matrix corresponding to the £1 transition "m is used 
in the nominal fit. To estimate the uncertainty, contri- 
butions from the quadrupole (M2) and other high order 
multipoles to the matrix are included in the fit, and 
the changes in the final results are treated as a systematic 
uncertainty. The uncertainties associated with the reso- 
nance parameters of intermediate states are estimated by 
varying their values by lo of their uncertainties quoted 
in the PDG [13. 'The uncertainty due to the barrier ra- 
dius when calculating Br(r) in Eq. (1] is estimated 
by alternative fits with r — 0.25 or 0.75 fm, respectively, 
where r — 0.6 fm is the nominal value. The uncertainty 
associated with the direct three-body decay y. > KK 


is estimated by alternative fits with other spin-parity hy- 
potheses, e.g. a 07 or 3^ non-resonant component in the 
KK or Kr systems. The largest changes in the signal 
yields are taken as systematic uncertainties. Assuming 
all the systematic errors are independent, the overall sys- 
tematic obtained by taking the quadrature sum of the 
individual values. 

In summary, the HSR. suppressed processes of 7-2 — 
K*(892)K and xc» — p(TT0)*m* are studied with the 
1/(3686) data collected by BESIII for the first time. 
The branching fractions of Xe — K*(892)*AT and 
Xe2 — K*(892) K? + c.c. are measured to be (1.5 + 
0.1 + 0.2) x 1074 and (1.3 + 0.2 + 0.2) x 1074, respec- 
tively, which are rather sizeable with ri to those of 
the HSR conserving decay Xe2 — VV [17 [19]. These 
branching fractions are at least one M. ol magnitude 
larger than the upper limit of the branching fraction of 
Xe2 — p(T10)*m* (1.1 x 1079). It is worth noting that 
this phenomenon is anticipated by the HSR violation 
mechanism proposed in Ref. Namely, the HSR vi- 
olation in x, + K*(892)K occurs via the intermediate 
meson loops due to the large U-spin symmetry breaking, 
while that in Xe3 > p(770)*7F is due to isospin symme- 
try breaking. Due to the large mass difference between 
s and u/d quarks, the U-spin symmetry is broken more 
severely in comparison with isospin symmetry. This re- 
sults in the larger decay branching for y,? + K* (892)K 
than that for yc — p(770)*x*. The results are cru- 
cial for further quantifying the HSR violation mecha- 
nisms and also provide deeper insights into the un- 
derlying strong interaction dynamics in the charmonium 
energy region. 
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